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Executive  Summary 


High-Frequency  (HF)  antenna  design  is  challenged  by  advanced  HF  communication 
systems,  such  as  HF-MIMO  and  HF-NVIS  for  vehicular  units,  and  the  new  vector¬ 
sensing  HF  direction-hnding  systems.  These  systems  require  small,  wideband,  and 
efficient  HF  antennas  with  specihed  radiation  patterns.  Regardless  of  the  tradeoffs 
that  are  made  to  get  to  a  specihc  HF  antenna,  its  far-held  pattern  is  the  hnal  gover¬ 
nance  of  system  performance. 

Measuring  far-held  HF  patterns  is  challenging — HF  wavelengths  range  from  10  m 
at  30  MHz  to  100  m  at  3  MHz.  Typical  pattern  ranges  measure  the  electric  held  on 
a  hemisphere  with  a  radius  10  wavelengths  from  the  antenna.  Consequently,  an  HF 
pattern  range  requires  a  1-km  radius  and  necessitates  other  approaches.  One  approach 
measures  the  near  held  on  a  standard  pattern  range  and  extends  the  resulting  near- 
held  model  to  a  far-held  estimate. 

Tapan  Sarkar  wrote  a  substantial  technical  paper  on  “near-to-far”  transformation 
using  spherical  expansions  [12].  An  implementation  of  Sarkar’s  near-to-far  transfor¬ 
mation  is  assessed  in  this  report.  Phase  1  of  this  report  assesses  this  implementation 
on  the  following  antenna  models: 

•  Analytic  antenna  models  where  the  near  held,  far  held,  and  the  spherical  ex¬ 
pansions  are  computable  in  closed  form. 

•  Mini-NEC  antenna  simulations  where  the  near  and  far  helds  are  computable 
and  numerical  quadratures  obtain  the  spherical  expansions. 

This  implementation  does  recover  the  far  held  from  the  near  held — provided  the 
spherical  expansion  includes  a  sufficient  number  of  terms.  Phase  2  applies  this  ex¬ 
pansion  to  pattern-range  measurements  and  extends  the  implementation  to  estimate 
the  number  of  terms  and  account  for  measurement  noise. 
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IX 


1  HF  Antenna  Measurements 


HF  antenna  design  is  challenged  by  advanced  HF  commnnication  systems  and  devel¬ 
opments  in  HF  direction  finding.  Fignre  1  shows  a  recent  HF-NVIS  antenna  developed 
for  gronnd-mobile  applications.  Challenges  for  any  compact  HF-NVIS  antenna  are 
the  trade-offs  between  bandwidth,  pattern  (npward  radiation),  and  efficiency  with 
the  necessary  adaptive  tnning.  This  antenna  and  its  inevitable  tradeoffs  should  be 
verified  by  measurements. 


Figure  1:  Compact  HF-NVIS  antenna  for  vehicle  and  ground  forces  (Reproduced  by 
permission  of  Military  Hi-Q  [9]). 


Figure  2  shows  co-located  HF  antennas  for  direction  finding.  Rather  than  an  HF 
array  with  the  half-wavelength  spacing,  this  co-located  array  offers  direction  finding 
using  polarization  rather  than  spatial  diversity.  One  challenge  for  HF  vector  sensors 
is  array  calibration.  Array  calibration  for  a  co-located  HF  array  is  equivalent  to 
measuring  the  array’s  response  to  an  electric  field  for  each  azimuth,  elevation,  and 
polarization. 


Figure  2:  Compact  HF  direction-finding  antenna  (Reproduced  by  permission  of  the 
Institution  of  Engineering  &  Technology  [11]). 
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Consequently,  both  the  HF-NVIS  antennas  and  HF  vector  sensors  require  far- 
held  measurements  to  establish  antenna  performance.  Figure  3  shows  the  Antenna 
Pattern  Range  located  at  the  Model  Range,  Space  and  Naval  Warefare  Systems  Center 
Pacihc  (SSC  Pacihc).  The  electric  held  produced  by  an  antenna  placed  in  the  Antenna 
Pattern  Range  is  measured  over  the  upper  hemisphere  at  a  radius  of  approximately 
80  ft  (24  m).  For  3-MHz  HF  antennas,  the  Antenna  Pattern  Range  measures  their 
near  helds  at  24  m  or  1/4  wavelength.  For  30-MHz  HF  antennas,  their  near  helds  are 
measured  at  2.4  wavelengths.  This  variation  in  the  near-held  measurements  raises 
the  following  question:  Is  the  near-to-far  transform  stable  with  respect  to  Antenna 
Pattern  Range  measurements  ranging  over  0.25  to  2.4  wavelengths? 


Figure  3:  Antenna  Pattern  Range  at  SSC  Pacihc. 

This  report  assesses  a  near-to-far  transformation  with  the  goal  of  mapping  the 
near-held  measurements  taken  on  a  grid  of  azimuth  and  elevation  angles  on  the  An¬ 
tenna  Pattern  Range  into  far-held  patterns.  Because  the  far-held  pattern  is  not 
directly  measurable,  this  Phase  1  discussion  relies  on  bounding  the  error  in  analytic 
models  before  undertaking  the  leap  to  the  Antenna  Pattern  Range  measurements. 
Accordingly,  Section  2  reviews  the  near-to-far  transformation  and  its  discrete  im¬ 
plementation  with  the  errors  arising  from  truncation  of  the  spherical  sums  and  the 
discrete  sampling  of  the  electric  held. 
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Section  3  reviews  Tapan  Sarkar’s  near-to-far  transformation  [12].  The  spherical 
expansions  and  the  special  functions  are  made  explicit  with  correct  normalizations, 
end  points  of  the  expansions,  and  three-term  recurrence  relations  are  supported  for 
computational  processing. 

Sections  4-11  apply  the  near-to-far  transform  to  selected  antenna  models.  Each 
section  undertakes  the  near-held  and  far-held  models  of  an  antenna,  an  “analysis” 
of  the  near  held  in  spherical  modes  using  only  discrete  samples  of  the  electric  held, 
and  a  comparison  the  resulting  synthetic  spherical  expansion  against  the  modeled 
antenna  in  both  the  near  and  far  helds.  These  examples  show  that  controlling  the 
truncation  error  is  more  critical  than  the  discretization  errors.  That  is,  sampling  the 
electric  held  at  5°  to  18°  increments  (discretization  error)  contributes  less  error  than 
omitting  higher-order  terms  in  the  spherical  expansion  (truncation  error). 

The  hnal  section  outlines  the  tasking  for  Phase  2  where  the  near-to-far  transfor¬ 
mation  will  be  applied  to  measurements  collected  on  the  Antenna  Pattern  Range. 
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2  Near-to-Far  Field  Transformations 


A  canonical  problem  in  antenna  engineering  is  estimating  the  inaccessible  far  field 
from  near-held  measnrements.  This  report  assesses  one  implementation  of  Tapan 
Sarkar’s  “near-to-far”  transformation  [12],  Figure  4  illustrates  this  near-to-far  imple¬ 
mentation. 


E(c/,0^,,d/) 


E(r,0,(l)) 


Analysis 


!■' . ! 


Synthesis 


Figure  4:  Computational  chart  for  the  near-to-far  transform. 


Discrete  near-held  samples  of  the  electric  held  are  collected  on  a  sphere  of  radius 
r  =  a.  The  “analysis”  box  uses  these  samples  to  compute  electric-held  coefficients 
{sm,n}  by  numerical  quadratures.  The  main  source  of  error  in  these  coefficients  {sm,n} 
arises  from  the  discrete  samples  because  the  numerical  quadratures  operate  on  two- 
dimensional  splines  of  the  electric  held  samples.  That  is,  the  coefficients  {sm,n}  are 
accurate  representations  of  an  electric  held  recovered  from  the  sampled  electric  held 
[6]. 

The  “synthesis”  box  computes  the  electric  held  at  any  radius  r  >  a  as  a  sum  of 
Sarkar’s  orthonormal  basis  functions  {ern,n}'- 

oo  n 

E(r,  e,(j))  =  0)  =:  </*)■  (1) 

n=0  m=0 

The  key  observation  is  that  the  expansion  holds  regardless  of  the  measurement  sphere: 

6(ai;  r,  9,  (/>)  =  6(02;  r,  9,  0) 


for  r  >  Oi,  02. 


Assessing  the  quality  of  this  expansion  focuses  not  only  on  computing  the  far-held 
pattern  using  a  truncated  sum 

N  n 

E(r,  6',  0)  ~  ^  ^  Sm^nGmA^,  9,  0)  =:  ©(o;  N;  r,  9,  0), 


n=0  m=0 
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but  a  truncated  sum  where  the  coefficients  of  the  expansion  have  accuracy  limited  by 
the  discrete  sampling: 


N  n 

E{r,e,(j))  ^  ^  ^Sm,n(A6',  A0)em,n(r,6',  0) 

71=0  m=0 

=  :  &{a;  A9,  Acj),  N;r,9,  (j))  (2) 

Thus,  the  expansion’s  accnracy  is  degraded  by  truncation  error  and  limited  by  the 
discretization  errors. 

Phase  1  of  this  report  assesses  Sarkar’s  expansion  on  several  antenna  models.  The 
Hertzian  dipole  is  the  baseline  computational  vehicle  for  Phase  1.  The  goal  is  to 
verify  the  numerical  computation  and  the  sensitivity  of  the  approximation 

E(r,  9,  0)  ^  6(a;  A9,  A0,  N;  r,  9,  0) 

on  this  analytic  dipole  model.  Once  this  ideal  dipole  is  analyzed,  other  antennas 
models  are  treated  in  the  same  format:  analysis,  synthesis,  and  the  far-held  pattern. 
Figure  5  summarizes  this  process  where  the  “ground  trnth”  of  the  antenna  model’s 
electric  held  is  compared  against  the  synthetic  held. 


Antenna 

Far  Field 

Model 

Near-Field 
Samples  ofE 


(liscretizcitum  errors  trwicutum  errors 


Figure  5:  Flow  chart  of  the  near-to-far  comparison  using  an  antenna  model. 
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Table  1:  Notation. 


Variable 

Description 

k 

Wavenumber 

A 

Wavelength  (m) 

/ 

Frequency  (Hz) 

u 

Frequency  (rad/sec) 

V 

Free  space  impedance  1207r  (ohms) 

e 

Permittivity  (F/m) 

r 

Radius  (m) 

e 

Zenith  (radians) 

4> 

Azimuth(radians) 

E 

Electric  held  vector  (V/m) 

;  = 

Dehnition 
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3  Sarkar’s  Near-to-Far  Transform 


Figure  6  illustrates  spherical  coordinates  and  the  decomposition  of  the  electric  field 
in  the  6  X(j)  plane.  Sarkar  expands  the  electric  field  using  spherical  Bessel  functions  of 
the  second  kind  hn  \kr),  associated  Legendre  functions  of  the  first  kind  P^{cos{6)) 
with  argument  cos(6'),  and  coefficients  amn,  Pmn,  'Jmn,  and  5mn  determined  by  the 
near-field  boundary  [12,  Equation  5]: 


Ee{r,  e,  0)  =  jcup  ^  ^  ^  ilran  sin(m0)  -  5ran  cos(m0)} 


n=0  m=0 


cxD  n  ^  Q 


sin(6*) 
dP^{cos{9)) 


n=0  m=0 


dd 


{a^n  cos(m0)  -f-  [3ran  sin(m0)}. 


and  [12,  Equation  6] 


°°  "■  d  P™- ( crm(  f)]] 

E^{r,e,(p)  =  jujfi  EE  hn\kr) — - —{imn  cos(m0)  5mn  sin(m0)} 


n=0  m=0 


oo  n  Q 


n=0  m=0 


dO 

mP^{cos{9)) 

sm{9) 


{a^n  sin(m0)  -  fd^n  cos(m0)}. 


Figure  6:  Spherical  coordinates  for  the  electric  field  showing  the  tangent  vectors  {E^ 
is  omitted). 
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The  coefficients  are  computed  from  the  near  held  at  radius  r  =  a  [12,  Equations 
14  and  27]: 


_ _ a _ (2n  +  l)(n  -  m)\  r 

kahn^'{ka)  +  hn\ka)  2n(n  +  l)(n  +  m)\  2ti  Jq  Jq 

Eg{a,  9,  0) — ^  sin(6*)  cos(m0)  —  mE^{a,  6,  0)P™(cos(6*))  sin(m0)  >  , 


[12,  Equations  14  and  Eq.  31]: 

a  (2n  +  l)(n-m)!  r 

kah^^\ka)  +  h^^\ka)Mn+l){n  +  m)\2^l 

I Ee{a,  9,  0)— sin(6')  sin(m0)  +  mE^{a,  9,  0)P™(cos(6'))  cos(m0)  1 , 


[12,  Equations  13  and  32]: 

^  1  (2n  +  l)(n-m)!  r 

jujfih^a\ka)  2n{n  +  l){n  +  m)\  27r  ^  Jq 

I Efi,{a,  9,  0)  sin(6')  cos(m0)  +  mEg{a,  9,  0)P™(cos(6'))  sin(m0)  1 , 


[12,  Equations  13  and  33]: 
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mn 


--  ^  {2n+l){n-m)\  7^^  r 

juj^ih^n\ka)  2n(n  +  l)(n  +  m)!  27r  Jq  ^  Jq 

E^{a,9,(j)) — — - ^sin(6')  sin(m0)  —  mEg{a,9,(j))P^{cos{9))  cos{m(j))  >  . 


The  twiddle  factor  is  the  standard  adjustment  for  working  with  the  sinusoidal 
basis: 


e 


m 


1  m  =  0 

2  m  >  0 


The  spherical  Bessel  functions  of  the  second  kind  are  dehned  in  Arfken  as  [3, 
Equation  11.150]: 

■=  (3) 

where  Hn  \r)  is  the  Hankel  function  of  the  second  kind.  Figure  7  plots  the  hrst  few 


spherical  Bessel  functions  [3,  Equation  11.150]: 


ho\r)  =  exp(-jr)- 

r 

hf\r)  =  exp(-jr)|-i  +  ^| 
hf{r)  =  exp(-jr)|-^-^  +  ^|. 


Figure  7:  Spherical  Bessel  functions  of  the  second  kind. 

The  associated  Legendre  functions  of  the  hrst  kind  P^{x)  are  [3,  Equation  12.84 
and  footnote]: 

Jm 

PT(x)  :=  (-ir(l  -  xT''^—P„(x), 

where  Pn{x)'s  are  the  Legendre  polynomials.  There  is  the  handy  observation  that 

P“(l)  =  Pn{x). 

A  stand-alone  form  adapted  from  [2,  Equation  4.102] 

n-ix)  :=  i^(l  -  xT'^—  {(x^  -  D”)  (4) 
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that  computes 


P^{x)  =  0  for  m  >  n  (5) 

Table  2  lists  the  initial  Legendre  polynomials  and  their  derivatives.  Figure  8  plots 
the  hrst  few  Legendre  polynomials. 


Table  2:  Legendre  polynomials  [3,  Table  12.1]. 


n 

Pn{x) 

PLiA 

0 

1 

0 

1 

X 

1 

2 

|(3a;^  -  1) 

3x 

3 

|(5a;^  —  3a;) 

i(15a;2  -3) 

4 

|(35a;^  -  30a;2  +  3) 

i(35a;3  -  15a;) 

Legendre  Polynomials 


Figure  8:  Legendre  polynomials. 

An  ambiguity  at  n  =  0  occurs  because  amn,  Pmn,  Imn  and  6mn  contain  a  “1/0” 
term  for  n  =  0.  If  these  coefficients  are  finite,  the  contribution  of  n  =  0  to  the 
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near-to-far  transformation  is  zero: 

juJiim-^^—^P^{cos{9)){'yrnn  sin(m0)  -  5mn  cos(m0)} 


n=0  m=0 


sin(6') 
r  (9r 


0  1  n 

EE-7);: 

n=0  m=0 


(iP™(cos(6*)) 


de 


{amn  cos(m0)  +  /3mn  sin(m0)} 


r  (9r 
=  0. 


1  <9  /  (2)n  \\  dPo{cos{6)) 

rPn\kr)) 


The  last  equality  follows  because  the  derivative  of  Pq{x)  =  Po(a^)  =  1  is  zero.  There¬ 
fore,  the  near-to-far  transformation  starts  the  double  sum  at  n  =  1. 


3.1  dP^  {cos{9))  /  dO 


Differentiating  the  generating  function  [3,  Equation  12.84]  computes  this  derivative 
in  terms  of  these  functions.  Recall  that  P^ix)  =  Pn{x)  denotes  the  Legendre  polyno¬ 
mials.  If  0  <  m  <  n, 


dPTi^] 


dx 


X 


=  -m - -i-lpil-x^p/^^Pnix) 

l-x^^  ’  ^  ’  dx^  ^  ’ 


_(_l)m+l(i  _  _  ^2^(m+l)/2 


^m+1 


■Pnix) 


=  -  (1  -  xY"Pp:‘*^(x). 

1  —  x^ 

In  terms  of  the  elevation  angle, 

dP^  (cos(6*) )  dP^  (x)  dx 

dO  dx  dO 

-m-^^PZ(x)  -  (1  -  w|  (-sta(9)) 

1  —  x^  J 

=  <  m^.^l[^lpP(cos(d))  +  -^P™+^(cos(0))|  sin(0) 

sm  (6^)  smlt^j  J 

=  mcot(0)Pr(cos(0))  +P™+^(cos(0)). 
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If  m  =  n,  the  (n  +  1)  derivative  of  Pn{x)  is  zero  using  =  0  from  Equation  4. 

Consequently, 

r/P™  frosts')') 

"  ^  ”  =  mcot(«)ir(™s(9))  +  C‘+'(cos(9))  (6) 

holds  for  all  0  <  m  <  n.  When  this  derivative  is  weighted  by  the  sine  function, 

sin(6')  =  mcos(6*)P™(cos(6*))  +  sin(6')P™’^^(cos(6')),  (7) 

du 

3.2  d/dr  (^rhn\kr)^ 

The  radial  term  is  computed  from  the  spherical  Hankel  function  of  the  second  kind 
using  the  three-term  recurrence  relation  [1,  Equation  10.1.12]: 

h(f)'(r)  =  hI2i(r)  - 

r 

Consequently,  the  derivative  of  the  spherical  Hankel  functions  can  be  computed  using 
these  functions  themselves: 

^  {rh//\kr))  =  h//\kr)  +  krh//'^' {kr) 

=  h//\kr)  +  kr  |/iI2i(/cr)  - 

=  krh^^/d^r)  —  nh//\kr) .  (8) 

The  interpretation  of  [12,  Equation  14]  is  computed  as  follows: 

^  ^  d'h//\kr))\^^^  =  kah^^l^ika)  -  nh//\ka). 
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4  Hertzian  Electric  Dipole 


The  Hertzian  dipole  is  the  baseline  to  assess  Sarkar’s  spherical  expansion.  Fignre  9 
illustrates  spherical  coordinates  around  a  short  antenna  carrying  a  z-directed  current 
Ia  in  free  space  on  a  wire  of  length 


Figure  9:  Hertzian  dipole  in  spherical  coordinates:  r  =  rr{6,(j)). 


The  next  two  hgures  present  the  nonzero  magnetic  and  electric  helds  of  the 
Hertzian  dipole  [7,  Section  17.2]: 


Er{r,6,  0) 

Ee{r,e,(j)) 


2 


Air 

iAAi  ,^2 


-k^  sin(6*)  exp{—jkr) 

9 

rj  cos(6*)  exp(— j/cr) 


j  1 
+ 


kr  {krY  j  ’ 
1  J 


=  V 


2n 

iA^ij  2 


{krY  {krY  j  ’ 


Air 


k'^  sin(6*)  exp(— jfcr)  + 


kr  {krY  {krY  j 


(9) 
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Hertian  Dipole  Er{r,0):  Al  =0.1  (meters);  /=300  MHz 


radius  (meters):  unit  current 


Figure  10:  \Er{r,9,  (f))\:  ^  =  1  (A),  A  =  1  (m),  Ai  =  A/10  (m). 


Hertian  Dipole  Ee{r^0):  Al  =0.1  (meters);  /=300  MHz 


1  23456789  10 


radius  (meters):  unit  current 

Figure  11:  \Eo{r,9,  (j))]:  M  =  1  (A),  A  =  1  (m),  Ai  =  A/10  (m). 
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The  analysis  and  synthesis  of  this  Hertzian  dipole  consists  of  sampling  the  helds 
near  the  antenna  (r  =  a)  and  recovering  the  far-held  patterns.  With 

1 


Erir,9,(j))  =  O 


[{kr) 


Eg{r,6,(j))  =  sm{6)  exp{—jkr)-^  +  0 


dTT 


kr 


[{kr) 


as  r  3>  a,  the  pattern  function  for  the  electrical  held  is  [7,  Section  17.4] 

\iA\M 


E{6, 0)  :=  lim  r  ||E(r,  9,  0)||  =  p- 


Att 


-ksm{9). 


(10) 


4.1  Hertzian  Dipole  Expansions 

This  section  uses  the  electric  held  of  the  Hertzian  dipole  to  assess  Sarkar’s  expansion 
in  two  stages.  For  the  hrst  stage,  Sarkar’s  electric-held  expansion  is  computed  from 
the  physical  model  of  the  dipole.  The  Hertzian  dipole  has  electric  helds  (Equation  9): 

E,(r)  :=  + 

Eg{r,9,(j))  =  Eo{r)sm{9) 

E^{r,9,(j))  =  0. 

Second,  the  coefficients  determined  by  this  physical  model  must  match  the  coefficients 
computed  by  quadratures  of  the  electric  held  measured  on  any  near-held  sphere. 

The  hrst  stage  will  verify  that  the  Hertzian  dipole’s  coefficients  are  all  zero,  with 
the  exception  of 

ttoi  =  -Vk—. — •  (11) 

An 

Comparing  the  physical  model  with  Sarkar’s  expansion  forces  the  expansion  in  func¬ 
tion  of  0  to  vanish  or  that  the  sums  collapse  to  m  =  0: 

Eg{r,9,(j))  =  Eo{r)sm{9) 

Pn  icos{9)) 


oo  n 


n=0  m=0 
oo  n 


sin(6') 
dPff{cos{9)) 


n=0  m=0 


d9 


{7mn  sin(m0)  -  cos(m0)} 

{amn  cos(m0)  -1-  Prnn  sin(m0)} 


^  ^  \^  dPn{cOs{9)) 


n=0 

oo 


n=0 
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where  the  P„’s  are  the  Legendre  polynomials.  More  precisely,  amn  =  0,  jSmn  =  0, 
7mn  =  0,  5mn  =  0  for  771  >  0.  Table  2  provides  the  derivatives  of  the  Legendre 
polynomials  to  substitute  into  Eg-. 

Ee{r,9,(j))  =  aoi-— (rh^^\kr)'^  Pl{cos{9)){- sm{e)) 

+ao2^^  (rh^2\kr)'^  P2{cos{9)){-sm{9))  +  ... 

=  (^rhf\kr)ym{9) 

1  O 

~«02-Tr  3cos(6*)  sin(0)  +  . .. . 

r  or  \  / 

Matching  the  sinusoidal  basis  in  9  gives  that  this  electric  held  component  expands  as 

Eo{r,9,(j))  =  -aoi--^  (rh^^\kr)^  sm{9). 

The  spherical  Bessel  term  expands  as 


r  dr 


rh^^\kr)  \  =  A;exp(— j/cr)  ^  f—  + 


3 


3 


kr  {krY  (kr)^ 
Substituting  into  the  electric  held  component  gives 

3 


(12) 


E0{r,9,(l))  =  —aoism{9)kexp{—jkr) 


3  1 


kr  {krY  (kr)^ 


Comparison  with  the  analytic  formula  for  Eq  in  Equation  9  forces  the  closed-form 
solution  given  in  Equation  11.  That  is,  starting  with  the  physical  model  of  Eg,  the 
Hertzian  dipole’s  coefficients  are  all  zero,  excepting  ctoi.  To  complete  this  hrst  stage, 
it  is  necessary  to  verify  that  this  expansion  explains  using  Sarkar’s  expansion 
(ttmn  =  0,  (3mn  =  0,  7mn  =  0,  5mn  =  0  for  m  >  0):  [12,  Equation  6] 

"■  9P™‘(co^(9)) 

E^{r,9,(j))  =  jufi  ZE  hn\kr) — - —{'Jmn  cos(m0)  -h  Smn  sin(m0)} 


n=0  m=0 


oo  n  o 


d9 

mP™(cos(6*)) 


n=0  m=0 

CXD 


sin(6*) 


{a^nn  sin(m0)  -  Pmn  cos(m0)} 


n=0 

=  0, 


dPn(cOs(9)) 
- d9 - 


provided  the  7mn’s  vanish  entirely. 
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The  second  stage  verifies  that  the  coefficients  in  the  expansion  can  be  computed 
using  only  the  electric  fields  measured  on  a  given  sphere.  Sarkar  computes  aoi  using 
Ee  and  measured  on  a  sphere  of  radius  a  whereas  the  expansion  of  the  electric  fields 
are  independent  of  any  spherical  measurements.  Consequently,  Sarkar’s  formulas 
must  cancel  the  measurements  on  the  sphere.  Equivalently,  Sarkar’s  formulas  must 
be  independent  of  the  radius: 


aoi  — 


(2n  +  l)(n  —  m)! 


'>271 


dfj)  /  d9 


kahn^\ka)  +  hn\ka)  2n(n  +  l)(n  +  m)!  27r  Jq  Jo 
E0{a,  9,  0) — — - ^sin(6')  cos(m0)  —  mE^{a,  9,  0)P™(cos(6'))  sin(m0) 


d9 


m=0 

n=l 


kah^J^' (ka)  +  h^J\ka)  4  27r 
a  3  r 

kah^J^Jka)  +  h^J\ka)  4  Jo 

3  a  i 


^^x2ttx  [  Eg{a,9,(j)) 


dPi{cos{9))  . 


d9 


sm{9)d9 


4  kahf^  (ka)  +  h^J\ka)  Jo 
Equation  12  computes  the  denominator: 


Eo{a,  9,  0)(—  sin(6*))  sm{9)d9 

TT 

Eg{a,  9)  sin(6')^(i6*. 


kahf^\ka)  +  hf\ka)  =  ^ 

On  the  sphere,  the  electric  field 

E,{a,0)  =  ,;^Un(f>)eM-3ka)  +  T  _ 

integrates  as 


Eg{a,  9)  sm.{9Yd9 


iAMk 

7]— - exp(-jA;a) 

All  a 

iAMk 

7]— - exp(-jA;a) 

All  a 


3 

{kaY 


3 

{kaY 


sm{9Yd9 
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Combining  the  numerator  and  denominator 


ttOl 


4  {ka)  +  h^i\ka)  Jo 

3  a 


E0{a,  9)  siniOYdO 

iA^i  k 


4exp(-jA;a)  + 


V- 


An  a 


X  exp(— j/ca) 


—Tjk 


Ato 


j 

{kaY 


Consequently,  Sarkar’s  denominator  term  cancels  the  effect  of  the  radius  a  in  the 
measured  electric  held  to  produce  an  agi  that  is  independent  of  the  measuring  sphere. 
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4.2  Hertzian  Dipole  Quadratures 

The  Hertzian  dipole  has  electric  helds  (Equation  9): 

Eg{r,e,(j))  =  Eo{r)sm{e), 

E^{r,e,(t))  =  0. 

Figure  12  plots  the  magnitude  of  Eg  on  the  sphere  a  =  3A  over  the  angular  grid 
sampled  at  18°  increments. 


Herlizan  Dipole 


Figure  12:  Hertzian  dipole  \Eg{a,9,(j))\;  m  =  1  amp,  A  =  1  meter,  Ai  =  A/10  meter; 
a  =  3A,A0  =  18°,  A0  =  18°. 

The  Hertzian  dipole  has  only  the  single  non-zero  coefficient  given  by  Equation  11: 

i^Ai 

aoi  =  —rjk - =  —18.8496, 

An 

where  the  numerical  value  is  computed  using  the  parameters  in  Figure  12.  Tables  3 
and  4  report  on  the  coefficients  amn  estimated  from  the  sampled  values  on  the  sphere 
a  =  3A.  The  real  part  of  these  amnS  is  accurate  to  5  places  whereas  the  imaginary 
part  is  accurate  to  6  places. 

Tables  3  and  4  compute  the  amn  from  the  electric  helds  sampled  on  a  sphere  of 
hxed  radius  a  =  3A.  Figure  13  displays  the  effect  of  varying  the  sampling  radius.  The 
accuracy  remains  at  5  places  for  A  <  a  <  5A.  This  stability  of  aoi  as  a  function  of 
the  sampling  radius  provides  a  partial  answer  to  the  question  of  Section  1.  The  near- 
to-far  transformation  is  stable  for  near-helds  of  the  ideal  dipole  where  the  sampling 
sphere  ranges  as  A  <  a  <  5A. 
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Table  3:  3fJ(Q!mn)  estimated  from  Eg  and  sampled  on  the  sphere  a  =  3A. 


n  =  1 

n  =  2 

n  =  3 

m 

=  0 

-18.84930 

0.00000 

-0.00000 

m 

=  1 

0.00000 

0.00000 

-0.00000 

m 

=  2 

-0.00000 

-0.00000 

m 

=  3 

-0.00000 

Table  4:  ^'(arnn)  estimated  from  Eg  and  E^  sampled  on  the  sphere  a  =  3A. 


n  =  1 

n  =  2 

n  =  2) 

m  = 

0 

-1.11  •  lO”^'^ 

-1.11  •  lO”^'^ 

-1.18-  10“® 

m  = 

1 

2.36-10-1^ 

-5.95  •  10-1^ 

-1.63-10-1® 

m  = 

2 

8.75  •  10-1® 

-5.65  •  10-1® 

m  = 

3 

-1.32  •  10-1® 

a  (meters) 

Figure  13:  Hertzian  dipole:  error  in  3fJ(aoi)  as  a  function  of  the  sampling  sphere. 


20 


4.3  Hertzian  Dipole  Synthesis 


Once  the  spherical  coefficients  are  computed  from  measurement  on  a  sphere  of  radius 
r  =  a,  Sarkar’s  expansions  compute  the  electric  held  at  distances  r  >  a.  This  section 
verihes  that  the  sampled  electric  held  can  be  recovered  on  the  measurement  sphere. 

Figure  14  reports  on  the  expansion  of  Eg{r^6^(f))  at  r  =  a  using  the  numerical 
coefficients  in  Tables  3  and  4.  Comparison  with  the  analytic  Eg{a,6,(j))  in  Figure  12 
shows  credible  agreement.  In  the  notation  of  Section  2,  the  synthetic  Eg{a,6,(j)) 
computed  using  Sarkar’s  expansion  is  a  credible  approximation  despite  the  truncation 
{N  =  3)  and  discretization  {A9  =  18°,  A0  =  18° ): 

Ee{a^  6,  0)  ~  &{a  =  3;  A9  =  18°,  A0  =  18°,  N  =  3;  a,  9,  0). 


Ee{r,  9,(p)  :r  =3  meters;  /  =300  MHz 


Figure  14:  Synthesis  of  Eg{a,  9,  0)  with  m  <  n  <  3. 
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Likewise,  Figure  15  reports  on  the  expansion  at  Erf){r,9,  cf))  at  r  =  a  using  the 
numerical  coefficients  in  Tables  3  and  4.  Because  E^[r,6,(j))  =  0,  the  absolute  error 
is  approximately  1.2  ■  10“^. 


E<j(r,  9,(p)  :r  =3  meters;  /  =300  MHz 


Figure  15:  Synthesis  of  E(i,{a,  0,  0)  with  m  <  n  <  3. 


4.4  Hertzian  Dipole  in  the  Far  Field 

The  last  step  in  Sarkar’s  near-to-far  transformation  is  to  take  the  coefficients  com¬ 
puted  in  the  near  held  and  synthesize  the  far-held  pattern.  Figure  16  shows  E0{r,  6,  (p) 
in  the  r-9  plane  out  to  r  =  lOOA.  Figure  17  shows  the  corresponding  image  of  the  syn¬ 
thesized  E0{r,  9,  (f))  determined  by  the  synthesis  code.  The  synthesis  code  implements 
Sarkar’s  expansion  given  the  coefficients  a,  /?,  7,  and  5.  This  code  was  tested  by 
inserting  the  exact  value  of  the  coefficients.  That  is,  all  coefficients  are  zeros  except¬ 
ing  aoi  =  —18.8496  set  by  Equation  11.  Comparison  of  both  hgures  shows  that  the 
magnitude  of  the  synthesized  electric  held  is  a  credible  approximation  to  the  analytic 
held. 
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Hertzian  Dipole;  /  =300  MHz 


Figure  16:  Exact  E0{r,9,(j)). 


Sarkar’s  Hertzian  Dipole:  /  =300  MHz 


Figure  17:  Synthetic  i76i(r,  0,  0). 
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The  far-field  pattern  function  for  the  electric  held  is  given  in  Equation  10  and 
repeated  here  for  convenience: 

F{9,(j))  :=  lim  r  ||E(r,6*,  0)||  =  — -ksm{6). 

r— >00  4:7T 

Figure  18  compares  this  pattern  function  against  two  approximations.  The  approxi¬ 
mation  labeled  “Exact”  is  r\EQ{r,9,  (j))]  evaluated  at  r  =  100  m.  The  approximation 
labeled  “Synthetic”  is  the  corresponding  computation  using  the  synthesis  code.  The 
approximations  are  relatively  close  to  the  true  pattern. 


Hertzian  Dipole  Patterns  at  /  =300  MHz 
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5  The  Small  Loop  Antenna 

The  Hertzian  dipole’s  computations  required  only  Eq  because  E^f,  is  zero.  A  loop 
antenna  tests  the  near-to-far  transformation  when  Efp  is  non-zero  and  Eq  is  zero.  The 
electric  held  for  a  small  loop  antenna  is  [4,  Equation  5-19] 

Er{r,9,(l))  =  0 

Ee{r,9,4))  =  0 

E^{r,9,(l>)  =  ,^!^A;3sin(0)exp(-jTr)|^--^|,  (13) 

where  AS*  =  Tirj  is  the  area  of  the  loop.  This  electric  held  expanded  and  synthesized 

to  assess  Sarkar’s  expansion. 


5.1  Expansions 


This  analytic  form  of  the  electric  held  determines  the  coefficients  of  Sarkar’s  expan¬ 
sion.  Sarkar’s  expansion  for  is  repeated  here  for  convenience:  [12,  Equation  6] 


E^{r,9,(l))  =  jA/i  EE  h<i\h 


n=0  m=0 


dP,["(cos(^)) 

d9 


{imn  cos(m0)  5mn  sin(m0)} 


-  EE)|  (^r ^  ( /cr ) )  {o^mn  sin(?Tl0)  /3mn  COs(?Tl0)}. 


sin(6') 


n=0  m=0 

The  symmetry  of  E^  with  respect  to  0  forces  the  sums  to  collapse  to  m  =  0: 


7-.  /  /I  ^dPn(cOs(9)) 

E^{r,9,(j))  =  jujfx^hl^\kr) - — — ^7o„. 

n=0 

More  precisely,  this  expansion  forces  =  0,  (3mn  =  0,  =  0,  and  5mn  =  0  for 

m  >  0.  Consequently,  the  expansion  for  Eq  must  take  the  form  [12,  Equation  5]: 


Ee{r,9,(f))  =  jufi 


n=0  m=0 


oo  n  ^  ^ 


PnicOs{9)) 
sm{9) 

dP^{cos{9)) 


{7mn  sin(m0)  -  Smn  cos(m0)} 


n=0  m=0 


d9 

O^On* 


{amn  cos(m0)  -1-  Prnn  sin(m0)} 


n=0 
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The  vanishing  Ee  =  ^  forces  amn  =  0  for  all  m,  n  >  0.  Returning  to  the  expansion 
of  Etj,,  the  product  of  the  radial  term  and  elevation  term  for  n  =  1  is 


=  exp(— jr) 


--  +  4  K-sin(0)) 


=  sin(6*)  exp(— jr)  - ^ 


/ 


Comparison  with  the  physical  model  of  E^  of  Equation  13  forces  all  7on’s  to  zero  with 
the  exception 

That  is,  the  physical  model  of  the  small  loop  antenna  forces 


jcn/i7oi  =  Tj- 


An 


For  convenience,  recall  that  [5,  Equation  9-14],  [5,  Equation  9-28],  [5,  Equation  10-33]: 


Co  = 


Solving  for  the  coefficient 


7  = 


k  = 


UJ 


Co 


7oi  =  - 


1 7]  iA^S 


=  -J 


j  fi  Att  00 

1  iA^Sk^ 


-3  Co 

-3 


y/JE  An  Co 
iAASP 


An  Co 
k^ 


Ja^S  ,  2 


An 

That  is,  the  small  loop’s  coefficients  are  all  zero,  excepting 


7oi  =  -3 


Aa^S  ,2 


An 


k\ 


(14) 


The  second  task  is  to  verify  that  the  coefficients  in  the  expansion  can  be  computed 
using  only  the  electric  helds  measured  on  a  given  sphere.  Sarkar  computes  701  using  Eq 
and  Efj,  measured  on  a  sphere  of  radius  a  whereas  the  expansion  of  the  electric  helds 
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are  independent  of  any  spherical  measurements.  Consequently,  Sarkar’s  formulas 
must  cancel  the  measurements  on  the  sphere.  Equivalently,  Sarkar’s  formulas  must 
be  independent  of  the  radius: 


7oi  = 


(2n +!)(„-»)!  ^ 


jujnh^n\ka)  2n(n  + l)(n  +  m)!  27r  Jo  Jo 
dP^{cos{9))  . 


E^{a,e,(j))- 
1 


d9 

3  1 


sin(6')  cos(m0)  +  mEe{a,  9,  4>)P^{cos{9))  sin(m0) 
dP^(cos(9))  . 


ju}jj,h^J\ka) 

4  27r 

1 

3 

j(jfih^i\ka) 

4/0 

1 

3  / 

d(f)  /  d9E^{a,9,(j))- 


d 


d9 


sin(6*) 


E^{a,  9,4))  cos(6')  \  sin(6')  d9 


jojfih^iJka)  ^  Jo 
1 


£"0(0,  9,  4)  sin(6*)^  d9 


jojjjihi^Jka)  47r 
1 


ka  {kaY  J  4  Jq 


sin(0)3  d9 


jojlih^J\ka)  '  ^TT 
rjk^  lA^S 


V-^^k^{.-hi\ka)) 


jUfJ,  4:71 
=  -jk: 


.,2'i'A^S 


477 


m=0 

n=l 


That  is,  7oi  determined  by  the  small  loop’s  physical  model  (Equation  14)  may  also  be 
computed  by  a  quadrature  that  integrates  the  electric  helds  measured  over  a  sphere. 
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5.2  Quadratures 

The  small  loop  antenna  has  electric  helds  E0{r,  6,(j))  =  0  and 

E^{r,e,ct>)  =  77^^fc3sin(0)exp(-j/cr)  |^  -  , 

where  AS*  =  vrr^  is  the  area  of  the  loop.  Fignre  19  plots  the  magnitnde  of  E^  on  the 
sphere  a  =  3X  over  the  angnlar  grid  sampled  at  18°  increments. 

Small  Loop  Antenna 


Fignre  19:  \E^{a,  9,  0)|  of  the  small  loop  antenna;  =  1  amp,  A  =  1  meter,  re  =  A/10 
meter;  a  =  3A,  A9  =  18°,  A0  =  18°. 

The  small  loop  antenna  has  only  the  single  non-zero  coefficient  given  by  Eqna- 
tion  14: 

7oi  =  =  -jO. 098696, 

dvr 

where  the  nnmerical  valne  is  compnted  using  the  parameters  in  Figure  19.  Tables  5 
and  6  report  the  coefficients  7^^  estimated  from  the  sampled  values  on  the  sphere 
a  =  3A.  The  real  part  of  these  'ymnS  is  accurate  to  9  places  whereas  the  imaginary 
part  is  accurate  to  5  places.  Figure  20  displays  the  effect  of  varying  the  sampling 
radius.  The  accuracy  remains  at  hve  decimal  places  for  A  <  a  <  5A.  This  stability 


Table  5:  'Slipfmn)  estimated  from  Eg  and  sampled  on  the  sphere  a  =  3\ 


n  =  1 

n  =  2 

n  =  3 

m 

=  0 

1.65  •  10-^'' 

4.36-  10-^« 

6.07-  lO-o 

m 

=  1 

1.34-10-20 

-3.66  •  10-12 

-4.40  •  10-20 

m 

=  2 

-9.34  •  10-20 

-3.04-10-1® 

m 

=  3 

6.20-10-21 

Table  6:  9(7mn)  estimated  from  Eg  and  E^  sampled  on  the  sphere  a  =  3X 


n  =  1 

n  =  2 

n  =  3 

m 

=  0 

-0.098695 

0.000000 

-0.000000 

m 

=  1 

-0.000000 

-0.000000 

0.000000 

m 

=  2 

-0.000000 

0.000000 

m 

=  3 

-0.000000 

of  7oi  as  a  fnnction  of  the  sampling  radins  provides  a  partial  answer  to  the  qnestion 
of  Section  1 — the  near-to-far  transformation  is  stable  for  the  near-helds  of  the  small 
loop  antenna  where  the  sampling  sphere  ranges  as  A  <  a  <  5A. 


Fignre  20:  Small  Loop  Antenna:  error  in  A(7oi)  as  a  fnnction  of  the  sampling  sphere. 
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5.3  Synthesis 


Synthesis  uses  the  coefficients  estimated  by  the  quadratures  to  reconstruct  the  helds  of 
the  small  loop  antenna.  Figure  21  plots  synthesized  from  the  estimated  coefficients 
using  only  the  relatively  coarse  sampling  of  the  E^f,  and  Eg  =  0  helds.  The  shape  and 
peak  values  are  in  agreement  with  Figure  19. 

E(f(r,  9,(f)):r  =3  meters;  /  =300  MHz 


Figure  21:  Small  Loop  Synthesis  \E^{a,9,(l))\]  =  1  amp,  A  =  1  meter,  ri  =  A/10 

meter;  a  =  3A,  A9  =  18°,  A0  =  18°. 
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5.4  Far-Field  Pattern 


This  section  compares  the  exact  far-held  patterns  of  the  small  loop  antenna  to  the 
synthesized  patterns.  Figure  22  plots  as  a  function  of  range  r  and  angle  9  from 
Equation  13.  Figure  23  plots  generated  using  the  synthesis  code  with  the  exact 
values  for  Sarkar’s  coefficients:  a,  /3,  7,  and  5  are  all  zero  excepting  701  =  —jO. 098696. 
These  slices  of  the  exact  and  synthetic  electric  held  verify  that  the  synthesis  is  accurate 
for  this  antenna. 


Small  Loop  Antenna  Es'-  f  =300  MHz 


10  20  30  ,  40  50  60,  70  80,  ,90  100 

r  (m):  iA=l  (A);  r^=0.1  (m) 

Figure  22:  Small  Loop  Antenna  \E^{r,6,(p)\. 


Sarkar’s  Small  Loop  Antenna  E4i  f  =.300  MHz 


Figure  23:  Synthetic  \E^{r,9,  (j))\. 
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The  pattern  function  for  the  small  loop  antenna  is  obtained  from  Equation  13: 
F{6,  (f))  =  lim  r||E(r,  6*,  0)||  =  sin(6*). 

r— >oo  47r 

Figure  24  compares  this  pattern  function  against  two  approximations.  The  approxi¬ 
mation  labeled  “Exact”  is  r|E6i(r,  0,  0)|  evaluated  at  r  =  100  m.  The  approximation 
labeled  “Synthetic”  is  the  corresponding  computation  using  the  synthesis  code.  The 
approximations  are  relatively  close  to  the  small  loop  pattern. 


Small  Loop  Patterns  at  /  =300  MHz 


Figure  24:  Far-held  pattern  comparison:  =  1  (A),  A  =  1  (m),  =  A/10  (m). 
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6  A  Synthetic  a  Field 


Define  a  synthetic  electric  field  with  an  =  1  as  the  only  non-zero  coefficient.  The 
angular  components  of  the  electric  held  are 


Eg{r,9,(l)) 

E^{r,e,(j)) 


dPi{cos{6)) 

Te 


cos(0) 


Pi^(cos(6*)) 

sin(6*) 


sin(0). 


Use  [3,  Table  12.3  with  footnote] 

pi(cos(0))  = -sin(0)  ^  =  _cos(f^). 

du 

Use  [3,  Equation  11.150]: 


hi  (r)  =  exp(— jr) 


Eg{r,9,(f))  =  -Eo(r)cos(6')cos(0) 
E^{r,9,(j))  =  +Eo{r)sm{(j)). 
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Figure  25  reports  the  electric  field  and  the  resulting  synthesis.  The  npper  panels 
are  sampled  using  =  73  and  Nq  =  43  sample  points  for  a  angular  grid  size 
of  A0  =  5°  and  A6  =  4.3°.  Only  these  sample  values  are  input  to  the  analysis 
code  to  estimate  Sarkar’s  coefficients.  The  lower  panels  show  the  resnlting  synthesis. 
Comparison  of  the  upper  panels  and  lower  panels  shows  that  the  magnitude  of  the 
synthesized  electric  is  a  credible  approximation  to  the  sampled  field. 


a^^=l:  /^300  MHz;  r=3  m;  N=3 
\Ee\  measured 


JKJ^ 

^nr 


0  100  200  300  0  100  200  300 


\E(j)\  measured 


\Ee\  synthesized 


J\JV 


0  100  200  300 


\E^\  synthesized 


0  100  200  300 


Figure  25:  Synthesis  of  the  an  =  1  field. 
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The  electric  field  was  synthesized  in  the  preceding  plots.  However,  there  is  the 
qnestion  if  the  electric  field  determined  by  the  an  is  represented  by  this  single  term. 

Question:  Is  an  the  only  non-zero  coefficient  of  this  electric  field? 

The  following  analysis  answers  this  qnestion  positively  and  leads  to  the  second 
qnestion  regarding  the  nnmerical  qnadrature. 

Question:  How  well  does  the  numerical  quadrature  estimate  an  compare  to  the 
vanishing  coefficients  of  this  electric  field? 

The  /dmn’s  vanish  because  of  the  orthogonality  between  cos((;/))  and  sin(0)  and  the 
basis  functions  sin(m0)  and  cos(m0),  respectively. 


Pmn 


(2n  -I-  l)(n  —  m)! 


kah^n^\ka)  +  h^n\ka)  2n{n  +  l){n  +  m)\  2n 


;mn 


+  bd-rnn}- 


Integrating  each  term  separately, 

p27T  piT 

b^-rnn  =  d(j)  d9  mE^{a,  9,  (j))P^{cos{9))  cos{m(j)) 

Jo  Jo 

d9  m{Eo{r)  sin(0)}P™(cos(6'))  cos(m0) 

/*27r  PTT 

=  mE^ir)  /  sin(0)  cos(m0)(i0  /  P^{cos{9))d9 
Jo  Jo 

=  0 


bo'^mn 


r  dP'^(cos(9)) 

d(j)  d9  Ee{a,9,<t>)  "  sin(^)  sin(m0) 


JO  JO  d9 

=  —EqIt)  /  d4>  d9  {cos{9)cos{4>)} — ^ sin(6*) sin(m0) 

Jo  Jo 

p27Z  P'1 

=  —E^ir)  /  cos(0)  sin(m0)(i0 

Jo  Jo 

=  0. 


d9 

dP^{cos{9))  . 


d9 


sin(6*)  cos{9)d9 


Thus,  the  (JmnS  must  vanish  in  a  credible  quadrature.  Table  7  reports  the  quadrature 
computations  (Jmn  from  Figure  25.  These  coefficients  vanish  nearly  to  the  machine  e. 


The  6mnS  vanish  mainly  because  the  orthogonality  of  the  sinusoids,  excepting  for 
m  =  1. 

^  1  {2n  +  l){n-m)\  e„  ^ 

Om.ri  —  /  -X  /  X .  ^9,mnj‘ 


ju^ih^n\ka)  2n{n  +  l){n  +  m)\  2n 
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Table  7:  |/9mn|  estimated  from  Eq  and  77^  sampled  on  the  sphere  a  =  3A. 


n  =  1  n  =  2  n  =  3 

m  =  0 

0  0 

m  =  1 

4.12-10-1^  9.51-10-1®  9.17-10-is 

m  =  2 

2.90  ■  10-1^  1.08  ■  10-is 

m  =  3 

2.47  ■  10-is 

d, 


(f)\mn 


r  dd  E^(a,  d,  Sin(d)  sin(m0) 


'*27r  ^TT 

d0  /  dd  {i7o(?")  sin(0)} 


d9 

dP™(cos(d))  . 


dd 


sin(d)  sin(m0) 


r*27r 


^o(^)  /  sin(0)  sin(m(?!))(i0 


dd 


^o(?^)7rx(?^  =  1) 


dPi(cos(d))^^ 


dd 


sin(d)dd 


by  parts 


Eo{r)Tix{m  =  1)  ^  Pm(cos(d))  sin(d)  T  —  /  cos(d)P^(cos(d))dd 


— Po(^")7J'x(^  =  1)  /  cos(d)P^(cos(d))dd. 


./o 


For  the  contribution  of  the  elevation  components, 
do-,mn 


27r  ^TT 

d0  /  dd  mP0(a,  d,  0)P™(cos(d))  cos(m0) 


/27r  ^TT 

dcj)  /  dd  {— Po(^")  cos(d)  cos(0)}P™(cos(d))  cos(m0) 


»27r 


=  mEolr)  /  cos(0)  cos(m0)d0  /  cos(d)P™(cos(d))dd 

Jo  Jo 

piT 

=  mEQ{r)TTxi'm  =  1)  /  cos(d)P,^(cos(d))dd 

Jo 

piT 

=  Eo{r)7rx{m  =  1)  /  cos(d)P,^(cos(d))dd. 

Jo 

Therefore,  d0;mn  —  d0;mn  =  0  and  the  dmn’s  must  also  vanish.  Table  8  reports  the 
quadrature  computations  from  Figure  25.  These  coefficients  vanish  nearly  to  the 
machine  e,  with  the  exception  of  612- 
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Table  8:  \5mn\  estimated  from  Eq  and  sampled  on  the  sphere  a  =  3A. 


n  =  1 

n  =  2 

n  =  3 

m  = 

0 

0 

0 

0 

m  = 

1 

1.82  ■  10-19 

3.89  ■  10-® 

9.76  ■  10-29 

m  = 

2 

9.01  ■  10-22 

9.73  ■  10-11 

m  = 

3 

6.28  ■  10-22 

Tables  9  and  10  report  the  magnitude  of  amn  and  7mn-  The  questions  raised  in 
this  section  are  both  answered  affirmatively  with  the  numerical  quadrature  and  the 
analysis. 


Table  9:  \amn\  estimated  from  Eq  and  E^  sampled  on  the  sphere  a  =  3A 


n  =  1  n  =  2  n  =  3 

m  =  0 

8.14 -lO-i^  2.81-10-9  2.22 -lO-i^ 

m  =  1 

1.00  2.83  ■  10-1^  1.42  ■  10-6 

m  =  2 

4.74  ■  10-^  3.55  ■  10-29 

m  =  3 

8.85  ■  10-9 

Table  10:  |7mn|  estimated  from  Eq  and  E^f,  sampled  on  the  sphere  a  =  3A 


n  =  1  n  =  2  n  =  3 

m  =  0 

7.22  ■  10-29  3.90  ■  10-29  2.32  ■  10-29 

m  =  1 

1.38  ■  10-29  4.91  ■  10-29  7.17  ■  IO-22 

m  =  2 

4.52  ■  10-22  1.56  ■  10-29 

m  =  3 

7.43  ■  10-22 
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7  A  Synthetic  S  Field 


Define  a  synthetic  electric  field  with  5i2  =  1  as  the  only  non-zero  coefficient.  The 
angular  components  of  this  electric  held  are 


Ee{r,e,(j))  =  cos{(t)) 

sin(6*) 

E^{r,d,(i))  =  jujfxhf\kr)^^‘^  sin(0). 

For  brevity,  introduce 


Eo{r)  ■=  juiih^2\kr)  =  ja;/rexp(-jA;r)  {  - 


+ 


J3 


kr  {krY  {krY 
using  [3,  Equation  11.150].  Use  [3,  Table  12.3  with  footnote] 

PlYosiO))  =  — 3cos(6*)  sin(6*) 


to  compute 


The  electric  helds  are 


dP^icosiO)) 

Te 


3  —  6cos(6*)^. 


Ed{r,  9,  0)  =  3Uo(r)  cos(6*)  cos(0) 
EYr,0,(l))  =  3Uo(r){l  -  2cos(6')^}sin(0). 
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Figure  26  reports  the  analysis  and  synthesis  of  this  electric  held  at  300  MHz.  The 
upper  panels  show  Eg{r,  6,  0)  and  E^{r,  6,  0)  sampled  on  an  angular  grid  at  r  =  3  m. 
The  angular  grid  uses  =  73  and  Ng  =  43  sample  points  for  an  angular  grid  size  of 
A(j)  =  5°  and  AO  =  4.3°. 


0 

0 


5i2=1:  /^300  MHz;  r=3  m;  N=3 


\Ee\  measured 


measured 


synthesized 


J\JK 

yr\r 


0  100  200  300 


Figure  26:  Synthesis  of  the  612  =  1  held. 

The  following  tables  report  the  magnitude  of  Sarkar’s  coefficients  estimated  from 
the  sampled  electric  helds.  All  coefficients  are  smaller  than  10“^  excepting  812  ~  1.02. 
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Table  11:  \amn\  estimated  from  Eq  and  sampled  on  the  sphere  a  =  3A  for  the 
5i2  =  1  antenna. 


n  =  1  n  =  2  n  =  3 

m  =  0 

6.21-10-1^  2.17-10-5  1.23-10-54 

m  =  1 

2.68  - 10-4  7.81  -  10-54  3.62  -  lO-^ 

m  =  2 

2.37-10-6  1.34  -  10-56 

m  =  3 

2.83  -  10-^ 

Table  12:  |/5mn|  estimated  from  Eq  and  E^f,  sampled  on  the  sphere  r  =  3A  for  the 
6i2  =  1  antenna. 


n  =  1 

n  =  2 

n  =  3 

m  = 

0 

0 

0 

0 

m  = 

1 

4.16  - 10-54 

2.25  -  10-55 

8.02  - 10-55 

m  = 

2 

2.49  -  10-55 

5.59  -  10-56 

m  = 

3 

4.97-10-55 

Table  13:  |7mn|  estimated  from  Eg  and  E^  sampled  on  the  sphere  r  =  3A  for  the 
6i2  =  1  antenna. 


n  =  1  n  =  2  n  =  3 

m  =  0 

3.85  - 10-55  2.49  -  10-5«  3.61  - 10-55 

m  =  1 

2.03  - 10-  55  7.40  -  10-  58  2.42  - 10-5® 

m  =  2 

7.97-10-59  5.15-10-58 

m  =  3 

2.22  - 10-59 

Table  14:  \5mn\  estimated  from  Eg  and  E^  sampled  on  the  sphere  r  =  3A  for  the 
(5i2  =  1  antenna. 


n  =  1 

n  =  2 

n  =  3 

m  = 

0 

0 

0 

0 

m  = 

1 

3.83  - 10-56 

1.02  - 10° 

1.22  - 10-55 

m  = 

2 

1.93  -  10-58 

4.13  -  10-9 

m  = 

3 

3.29  - 10-58 
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8  Horizontal  Electric  Dipole 


Figure  27  shows  an  electric  dipole  of  length  Ai  centered  at  the  origin  and  oriented 
along  the  x-axis. 


Figure  27:  Horizontal  Dipole. 

If  the  dipole  is  carrying  current  iA  amps,  the  near  field  is 

Eeir.O^cf))  =  i7i(r)  cos(6')  sin(26*)  sin(0)^  +  cos(0)| 

+F'o(^")  sin(6')  cos(6')^  sin(0)^  (15) 

E^{r,6,4>)  =  — i7i(r)  cos(6*)  sin(0)(cos(6*)  —  cos(0)  sin(6*)) 

+Eo{r)  sin(6')  sin(0)(sin(6')  +  cos(0)  cos(6')),  (16) 

where 

:=  +  (17) 

E,(r)  :=  (18) 


41 


8.1  Derivation 


One  way  to  obtain  the  near  field  of  a  horizontal  dipole  is  to  rotate  the  electric  field 
of  the  Hertzian  dipole.  The  Hertzian  dipole  has  electric  fields  [7,  Section  17.3]: 

EeiV;r,e,(j))  =  Eoir)sm{e) 

E^{V]r,e,(j))  =  0 
Er(y;r,9,(l))  =  Ei{r)cos{9), 

where  the  “H”  makes  explicit  that  these  fields  reference  the  vertical  dipole.  The 
right-handed  matrix  R  that  rotates  z  into  x  while  leaving  y  fixed  is 

■  0  0  1' 

R=  0  10. 

-10  0 

If  E(i7;  r,  9,  0)  denotes  the  horizontal  dipole  field, 

E(77;r,0,0)  =  RE{V;r,9,  <j)) 

=  R^^Er{V-,r,9,(j))r  +  Eg{V-,r,9,(j))d^ 

=  Er(y-r,9,(l>)R^+Eg(y-,r,9,<P)Rd 
=  Er{V;r,9,(l))  |(i?r,r)r -h  {Rr,6)6  +  {Rr,<p)<p^ 

+Eg{V-,r,9,(j))  |(i?0,r)r-h  {R6,6)6  +  {R6,(f))(f)^ 

=  r  ^^Er{V-,r,9,(j)){Rr,r)  +  Eg{V-,r,9,  (j)){Rd,r)^ 

+e  r,  9,  6)  +  Eg{V-,  r,  9,  <i>){Re,  6) } 

+$  r,  9,  <t>){Rx,  0)  +  r,  9,  <i>){Re,  0) }  . 

The  action  of  the  rotation  matrix  R  on  the  radial  vector 

0  0  1  sin(6')  cos(0) 

Rx  =  0  10  sin(6')  sin(0) 

—  10  0  cos(6*) 

cos(6*) 

=  sin(6*)  sin(0) 

—  sin(6')  cos(0) 

=  {Rx^y)t  +  {Rx,6)6  +  {Rx,4))4> 
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has  coefficients 


{Rr,r)  =  sin(6')^  sin(0)^ 

{Rr,6)  =  -  sin(26')  sin((;/))^  +  cos(0) 

{Kr,(f))  =  —  sin(0)(cos(6')  —  cos(0)  sin(6')). 

The  action  of  the  rotation  matrix  R  on  the  elevation  vector 

0  0  1  cos(6')  cos(0) 

RO  =  0  10  cos(6')  sin((;/)) 

—  10  0  —  sin(6') 

—  sin(6*) 

=  cos(6*)  sin(0) 

—  cos(6')  cos(0) 

=  {Re,?)?  +  {Re,  6)6  +  {Re, 

has  coefficients 

{Re,?)  =  -  sin(26*)  sin(0)^  —  cos(0) 

{Re,  6)  =  cos(6*)^  sin(0)^ 

{Re,  <p)  =  sin(0)(sin(6*)  +  cos(0)  cos(6')). 

The  elevation  component  of  the  horizonal  dipole  is 

Ee{H-,r,e,4>)  =  Er{V-,r,e,<t>){R?,e)  +  Ee{V-,r,e,<t>){Re,e) 

=  Er{V;r,6,  0)  sin(26')  sin(0)^  +  cos(0)| 

+Eg{V ;  r,  9,  0)  cos(9)^  sin(0)^ 

=  i?i(r)  cos(6')  sin(26*)  sin(0)^  +  cos(0)| 

+Eo(r)  sin(6*)  cos(9)^  sin(0)^. 

The  azimuth  component  of  the  horizonal  dipole  is 

E^{H-,r,e,<t>)  =  Er{V-,r,e,<t>){R?,^)+Ee{V-,r,e,4>){Re,^) 
=  —Er{V-,r,9,(i))?m{(j)){cos{9)  —  cos(0)  sin(0)) 
+Eg{y-,r,9,(j))  sin(0)(sin(6*)  +  cos(0)  cos(6*)) 
=  —Ei{r)  cos(6')  sin((;/))(cos(6*)  —  cos(0)  sin(6')) 
+EQ{r)  sin(6*)  sin(0)(sin(6*)  +  cos(0)  cos(6*)). 
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8.2  Analysis 


This  section  verifies  that  the  and  the  7mn’s  of  Sarkar’s  expansion  vanish.  The 

/Smn’s  vanish  because  of  the  orthogonality  between  cos(0)  and  sin((;/))  and  the  basis 
functions  sin(m0)  and  cos(m0),  respectively, 


/3mn 


(2n  +  l)(n  —  m)! 


kah^n^'ika)  +  h^n\ka)  2n{n  +  l){n  +  m)\  2tt 


T  ^0;mn}  • 


Integrating  each  term  separately, 

p27T  piT 

b^-mn=  /  d(j)  d9  mE^{a,9,(l))P^{cos{6))cos{m(j)) 

Jo  Jo 

r27T  riT 

=  d(j)  d9  m{—Ei{r)  cos{9)  sm{(j)){cos{9)  —  cos{(l))  sm{9)) 

Jo  Jo 

+EQ{r)  sin(6*)  sin(0)(sin(6*)  +  cos(0)  cos(0))}P™(cos(6*))  cos(m0) 

/*27r  riT 

=  —mEi{a)  /  sin(0)  cos(m0)(i0  /  cos{9YP^{cos{9))d9 

Jo  Jo 

p27T  piT 

+mEi{a)  /  sin(0)  cos(0)  cos(m0)(i0  /  cos(6*)  sin(6*)P™(cos(6'))(i6' 

Jo  Jo 

n27T  PTT 

+mPo(a)  /  sin(0)  cos(m0)(i0  /  sm{9yP^{cos{9))d9 

Jo  Jo 

r27T  piT 

+mEo{a)  /  sin(0)  cos(0)  cos(m0)(i0  /  sin(6')  cos(6')P™(cos(6*))(i6* 
Jo  Jo 

=  0 


because  the  cos(m0)’s  are  even  functions  while  sin((;/))  and  sin((;/))  cos^cj))  are  odd  func¬ 
tions. 


b 


6\mn 


dP^ico^iO\\ 

d9  Ee (a,  9,  0)  "  ^  sin(0)  sin(m0) 

I  du 

r  .  ,.,rf^r(cos(0)) 

cos(0j  sin(m0)a0  /  cos(c^)  sinlcyj - — - 

1  ./n  d9 

(iP™(cos(6*)) 


r*27r 


cos((?!))  sin(m0)(i0  /  sin(6*)  cos(6*) 


d9 

d9 


Thus,  the  PmnS  must  vanish  in  a  credible  quadrature  (see  Table  16). 
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The  7mn’s  vanish  because  of  the  orthogonality  between  cos(0)  and  sin((;/))  and  the 
basis  functions  sin(m0)  and  cos(m0),  respectively. 

1  (2n  +  l)(n  —  m)! 


Tmn  — 


jujfihia\ka)  2n(n  +  l)(n  +  m)! 
Integrating  each  term  separately, 

dP^{cos{9))  . 


{C0  ;mn  T  CQ-rnn\- 


/2n  pTV 

djcj)  I  do  E (p{(i ^  0 ^  (j)^ ~ 


dO 


sin(6*)  cos{m(j)) 


r27T  piT 

/  dcj)  dO  {—Ei{a)cos{9)sm{(l)){cos{9)  —  cos{(l))sm{9)) 

Jo  Jo 

dP™-(coH(9)) 

+£'o(a)  sin(6*)  sin(0)(sin(6*)  +  cos(0)  cos(6'))} — ^  sin(6*)  cos(m0) 


d9 


r*27r 


=  —Ei{a)  /  sm{(p)cos{m(p)d(p  /  sin(6*)  cos((9) 


2dP™(cos(0)) 


d9 


d9 


r-27r 


+Pi(a)  /  sin((;/))  cos(0)  cos(m0)(i0  /  cos(6')  sin(6') 


2dP™(cos(0)) 


d9 


d9 


'■271 


+Po(o)  /  cos{m4>)d4>  /  sin(6') 


3dP™(cos(0)) 


'0 


'0 


d9 


+Po(fl)  [  sin((;/))  cos(0)  cos(m0)(i0  f  sin(6')^  cos(6') — 


>0 


>0 


d9 


=  0 


because  the  cos(m0)’s  are  even  functions  while  sin((;/))  and  sin(0)  cos^cj))  are  odd  func¬ 
tions. 


^27r  piT 

d(j)  /  d9  mEo{a,9,(l))P^{cos{9))sm{m(l)) 


'*277  /*7r 


=  /  (i0  /  (i6' m{Pi(a)  cos(6'){sin(26*)  sin((;/))^/2 -|- cos((;/))} 

Jo  Jo 

-|-Po(a)  sin(6')  cos(6*)^  sin(0)^}P™(cos(6'))  sin(m0) 

=  -—Pi (a)  /  sin(0)^ sin(m0)(i0  /  cos(6') sin(26')P™(cos(6*))(i6' 

2  Jo  Jo 

r27T  piT 

+mEi{a)  /  cos(0)  sin(m0)(i0  /  cos{9)P^{cos{9))d9 

Jo  Jo 

p27T  piT 

+mEo{a)  /  sin(0)^  sin(m0)(i0  /  sin(6*)  cos(0)^P™(cos(6*))(i6* 

Jo  Jo 

=  0 

because  the  sin(m0)’s  are  odd  functions  while  cos(0)  are  sin(0)^  are  even  functions 
over  [0,  27r].  Thus,  the  7mn’s  must  vanish  in  a  credible  quadrature  (see  Table  17). 


45 


Table  15:  \amn\  estimated  from  Eg  and  sampled  on  the  sphere  a  =  3A  for  the 
Horizontal  Dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

1.90 

5.31  ■  10-® 

1.27 

1.21  ■  10“® 

1.28  ■  10-® 

m  = 

1 

9.48 

5.01  ■  10-^ 

1.41  ■  10“® 

2.48  ■  10“® 

1.88  ■  10-^ 

m  = 

2 

1.03  ■  10“^ 

0.212 

6.82  ■  10-9 

9.84  ■  10-® 

m  = 

3 

1.19  ■  10“® 

7.48  ■  10“® 

4.32  ■  10-19 

m  = 

4 

1.32  ■  10-9 

8.51  ■  10-9 

m  = 

5 

1.36  ■  10-19 

Table  16:  \Pmn\  estimated  from  Eg  and  E^  sampled  on  the  sphere  r  =  3A  for  the 
Horizontal  Dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

0 

0 

0 

0 

0 

m  = 

1 

8.31  ■  10-1® 

4.46  ■  10-1^ 

3.61  ■  10-1^ 

7.05  ■  10-1® 

5.58 

o 

1 

00 

m  = 

2 

1.74  ■  10-1® 

3.12  ■  10-1^ 

3.14  ■  10-1® 

1.40 

O 

1 

m  = 

3 

6.61  ■  10-1^ 

9.67-10-1® 

1.58 

.  10-18 

m  = 

4 

3.33  ■  10-1^ 

2.79 

■  10-19 

m  = 

5 

3.44 

.  10-18 

Table  17:  |7mn|  estimated  from  Eg  and  E^  sampled  on  the  sphere  r  =  3A  for  the 
Horizontal  Dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

1.36  ■  10-19 

1.35  ■  10-19 

1.53  -  10-19 

7.26 

- 10-19 

1.68 

-  10-19 

m  = 

1 

9.59  ■  10-19 

1.07-10-1® 

2.02  -  10-19 

8.53 

- 10-21 

2.78 

-  10-29 

m  = 

2 

3.85  - 10-29 

8.87  -  10-29 

6.03 

-  10-29 

2.41 

-  10-29 

m  = 

3 

1.89  -  10-21 

2.42 

-  10-29 

1.20 

-  10-29 

m  = 

4 

7.81 

-  10-21 

6.08 

-  10-21 

m  = 

5 

1.56 

-  10-21 
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Table  18:  \5mn\  estimated  from  Eq  and  sampled  on  the  sphere  r  =  3A  for  the 
Horizontal  Dipole. 
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8.3  Synthesis 


The  following  figures  show  the  synthesis  for  this  antenna.  The  electric  helds  are 
sampled  on  an  angular  grid  at  a  radial  distance  of  3  m.  This  angular  grid  is  determined 
by  sampling  N(j,  =  73  and  Ne  =  43  for  a  angular  grid  size  of  A0  =  5°  and  A6*  =  4.3°. 
Figure  28  is  a  3-D  view  of  the  magnitude  of  the  electric  helds.  The  upper  panels  show 
Eg{r,9,(j))  and  E^{r,9,(j))  sampled  on  an  angular  grid.  The  lower  panels  show  the 
corresponding  synthesis.  Figure  29  displays  the  same  information  in  two-dimensional 
images. 

Horizontal  Dipole:  f=300  MHz;  t^3  m 


Figure  28:  Horizontal  dipole  at  300  MHz  in  3-D;  m  =  1  amp,  i  =  A/10;  upper  panels 
are  the  simulated  helds,  lower  panels  are  the  synthesized  helds. 
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Horizontal  Dipole:  f=300  MHz;  r=3  m;  A^=5 


\E0\  measured 


\E(j\  measured 


150 


0 


100  200  300 


\Ee\  synthesized 


\E(j\  synthesized 


0  100  200  300 


Figure  29:  Horizontal  dipole  images  at  300  MHz;  m  =  1  amp,  i  =  A/10;  upper  panels 
are  the  simulated  fields,  lower  panels  are  the  synthesized  fields. 
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8.4  Far  Field 


From  Equations  17  and  18,  the  radial  terms  behave  like 

Eo{r)  =  kexp{—jkr)-  +  0[r~‘^] 

Ei(r)  =  0[r-^]. 

Substitution  into  Equations  15  and  16  gives  the  behavior  of  the  electric  held: 

Ee{r,  9,  0)  =  Eo{r)  sin(6*)  cos(6')^  sin(0)^  +  0[r~‘^] 

E^{r,6,4>)  =  Eo(r)  sin(6*)  sin(0)(sin(6*)  +  cos(0)  cos(6*))  + 

The  magnitude  squared  of  the  far  held  has  the  expansion 

||E(r,0,0)f  =  |E,(r,0,0)P  +  |E^(r,0,0)P  +  |E,(r,0,0)|2 
=  |E'o(r)p  sin(6*)^  cos(6*)"‘ sin(0)"^ 

+  |E'o(r)|^  sin(6*)^  sin(0)^(sin(6')  +  cos(0)  cos(6'))^  +  0[{kr)~^] 

=  |Eo(r)p  sin(6')^  sin(0)^  {cos(6*)"‘  sin(0)^  +  (sin(6')  +  cos(0)  cos(6'))^} 
+0[{kr)-^]. 

The  pattern  function  of  the  electric  held  is 
E{e,(j))  =  lim  r||E(r, 0,0)11 

r— >oo 

=  lim  r|Eo(r)|  sin(0)  |  sin(0)| 

r— >CXD 

X  a/ cos(0)"^  sin(0)2  +  (sin(0)  +  cos(0)  cos(0))2 

\iA\Aik  .  . 

=  V — 1 - sm(0)  sm(0) 

4:71 

X  a/ cos(0)"^  sin(0)2  +  (sin(0)  +  cos(0)  cos(0))2. 

Figure  30  plots  this  far-held  pattern.  Figure  31  is  the  approximation  to  the  pattern 

E{9, 0)  ^  ry/|E,(r,0,0)|2  +  |E^(r,0,0)|2. 

at  r  =  100  m.  Figure  32  plots  the  corresponding  approximation  to  the  pattern 
function  using  the  electric  held  synthesized  by  Sarkar’s  coefficients.  The  pattern 
functions  are  comparable  in  all  three  hgures. 
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Figure  30:  Pattern  function  for  the  horizontal  dipole:  M  =  1  amp,  £  =  A/10  m. 


(m 


Figure  31:  Pattern  function  approximation  at  r  =  100  m. 
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Figure  32:  Synthetic  pattern  function  approximation  at  r  =  100  m. 
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9  Tilted  Dipole 

Figure  33  is  a  schematic  of  a  tilted  dipole  raised  over  an  infinite  ground  plane. 


Figure  33:  Dipole  tilted  in  the  x  x  2;-plane  and  lifted  over  a  perfect  ground. 

This  dipole  is  simulated  in  Mini-NEC  with  the  following  specihcations: 

•  Length  £  =  8  ft 

•  Tilted  from  vertical  in  the  x  x  2;-plane  by  Ot  =  45° 

•  Lifted  over  the  ground  by  he  =  6  ft 

Table  19  lists  the  Mini-NEC  description  for  completeness,  where  the  dipole  is  driven 
at  29  MHz.  Figure  34  reports  the  resulting  electric  helds  measured  on  a  hemisphere 
with  a  radius  of  80  ft  or  24  m. 
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Table  19:  MiniNEC  description. 


D : \antenna\NearToFarField\TiltedDipole\TiltedDipole01_80  01-22-2013  09:51:19 

Tilted  Dipole 

45  degrees  elevation  in  the  x-z  plane 
8  long 

6  ft  above  the  ground 


GEOMETRY 

Dimensions  in  feet 
Environment:  perfect  ground 


wire  caps  X  Y 

1  none  2.8284  0 

-2.8284  0 

Number  of  wires  =  1 


current  nodes  =  10 


Z  radius  segs 

8.8284  .001  11 

3.1716 


minimum 

Individual  wires 

wire 

value 

segment  length 

1 

.727266 

segment/radius  ratio 

1 

727.266 

radius 

1 

l.E-03 

maximum 
wire  value 

1  .727266 

1  727.266 

1  l.E-03 


ELECTRICAL  DESCRIPTION 
Frequencies  (MHz) 
frequency 

no.  lowest  step 

1  29.  0 


no .  of 
steps 
1 


segment  length  (wavelengths) 
minimum  maximum 

.0214425  .0214425 


Sources 

source  node  sector  magnitude  phase 

15  11.  0 


type 

voltage 
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e  (deg)  e  (deg) 


0 


Ee\‘^+ 


0  100  200  300 


r  E|f  at  80  (feet) 


Ee 


0  100  200  300 

(j)  (deg):  Pj  =1.05e-05  Watts 


(j)  (deg):  29  MHz 


Figure  34:  MiniNEC  simulation  of  the  tilted  dipole  at  29  MHz;  electric  fields  are 
measured  at  80  ft. 


55 


9.1  Synthesis 


Synthesis  requires  extending  the  electric  helds  from  the  upper  hemisphere  to  the 
entire  sphere.  If  Eg{9,(j))  and  denote  the  simulated  electric  helds  on  the 

upper  hemisphere,  make  the  even  extension 

Ee{e,(j>)  = 

0  <  0  <  7r/2 

7r/2  <  9  <  TT 

and  the  odd  extension 

E^{9,(p)  =  1 

Et(e,4>) 

0  <  0  <  7r/2 

7r/2  <  9  <  71 

Figures  35  and  36  report  the  magnitude  of  the  sampled  fields.  The  mesh  plot  is  used 
to  emphasize  that  fields  are  sampled  on  the  upper  hemisphere  with  =  73  samples 
in  0  <  0  <  360°  using  the  azimuth  increment  =  5°  and  Ng  =  19  samples  in 
0  <  0  <  90°  using  the  zenith  increment  =  5°.  Figure  36  verihes  that  vanishes 
at  the  ground  plane. 


Tilted  Dipole 


Figure  35:  MiniNEC  simulation  of  the  tilted  dipole  at  29  MHz;  electric  helds  are 
measured  at  80  ft. 

Tables  20,  21,  22,  and  23  report  the  coefficients  estimated  from  the  sampled  helds. 
The  relatively  small  size  of  suggests  that  these  coefficients  are  zero.  Likewise, 
the  Pmn  and  5mn  are  relatively  small  with  respect  to  amn-  Consequently,  when  the 
synthesis  undertaken,  the  a  terms  will  dominate. 
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Tilted  Dipole 


Figure  36:  MiniNEC  simulation  of  the  tilted  dipole  at  29  MHz;  electric  fields  are 
measured  at  80  ft. 


Table  20:  \amn\  estimated  from  Eq  and  sampled  on  the  sphere  at  80  ft  for  the 
tilted  dipole. 


n  = 1  n  =  2  n  =  3 

m  =  0 

3.63  ■  10-^  1.09  ■  10“^  2.72  ■  10“^ 

m  =  1 

9.48  ■  10“^  6.35  ■  10-3  1.05  ■  10“^ 

m  =  2 

5.91 -10-®  3.30-10-3 

m  =  3 

8.49  ■  10-^ 

Table  21:  \(3mn\  estimated  from  Eg  and  E^  sampled  on  the  sphere  at  80  ft  for  the 
tilted  dipole. 


n  = 1  n  =  2  n  =  3 

m  =  0 

0  0  0 

m  =  1 

1.32  ■  10-3  1.95  ■  10-^  5.98-10-^ 

m  =  2 

1.34  ■  10-^  4.23  ■  10-3 

m  =  3 

4.31  ■  10-3 
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Table  22:  |7mn|  estimated  from  E0  and  sampled  on  the  sphere  at  80  ft  for  the 
tilted  dipole. 


n  =  1  n  =  2  n  =  3 

m  =  0 

4.04  ■  10-11  4.06  ■  10-11  1.80  ■  lO-n 

m  =  1 

5.80  ■  10-*^  1.93  ■  10-*^  6.25  ■  10-1° 

m  =  2 

4.00  ■  10-1°  1.65  ■  10-1° 

m  =  3 

5.45  ■  10-11 

Table  23:  \5mn\  estimated  from  Eq  and  E^f,  sampled  on  the  sphere  at  80  ft  for  the 
tilted  dipole. 


n  = 1  n  =  2  n  =  3 

m  =  0 

0  0  0 

m  =  1 

4.95  ■  10-^  6.20  ■  10-^  6.46  ■  10-^ 

m  =  2 

3.04  ■  10-°  3.92  ■  10-° 

m  =  3 

2.43  ■  10-° 
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Figure  37  compares  the  sampled  and  synthesized  helds  at  29  MHz  on  a  upper 
hemisphere  with  an  80-ft  radius.  The  upper  panels  are  the  sampled  helds  computed 
from  MiniNEC.  The  lower  panels  are  the  synthesized  helds  computed  from  Sarkar’s 
expansion  with  only  N  =  3  terms  and  the  coefficients  amn,  Pmn,  1mm  ^mm  estimated 
from  the  sampled  helds  in  Figures  35  and  36.  These  images  use  the  same  linear 
color  scale  with  the  magnitude  of  the  helds  plotted  in  V/m.  Comparison  between 
the  sampled  helds  from  MiniNEC  and  the  synthesized  held  of  Sarkar  shows  that 
diherences  in  the  magnitudes  are  relatively  small. 
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Figure  37:  Tilted  dipole  synthesis  at  29  MHz  and  the  80-ft  radius.  The  upper  panels 
are  the  sampled  helds  computed  from  MiniNEC.  The  lower  panels  are  the  synthesized 
helds. 
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9.2  Far  Field 


Figure  38  compares  the  far  field  computed  by  MiniNEC  to  Sarkar’s  expansion.  The 
upper  panels  show  the  magnitude  of  the  electric  helds  (V/m)  computed  by  MiniNEC 
at  the  500-ft  radius.  The  sampling  still  uses  the  5°  increments.  The  lower  panels 
show  the  magnitude  of  the  electric  helds  computed  by  Sarkar’s  expansion.  That  is, 
the  electric  helds  were  sampled  at  an  80-ft  radius,  Sarkar’s  coefficients  were  estimated 
from  these  samples,  and  these  coefficients  compute  the  held  at  the  500-ft  radius.  The 
magnitude  of  these  plots  show  good  agreement. 
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Figure  38:  Comparing  far-held  patterns  of  the  tilted  dipole  at  29  MHz.  (upper 
panels)  sampled  helds  computed  from  MiniNEC  at  the  500-ft  radius,  (lower  panels) 
synthesized  helds  computed  from  the  sampled  helds  at  the  80-ft  radius  and  extended 
to  the  500-ft  radius. 
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10  Shifted  Tilted  Dipole 


Figure  33  is  a  schematic  of  the  tilted  dipole  from  the  preceeding  section  that  is  shifted 
along  the  x-axis. 


A 


Figure  39:  Dipole  shifted  along  the  x-axis,  tilted  in  the  x  x  2;-plane  and  lifted  over  a 
perfect  ground. 

This  dipole  is  simulated  in  MiniNEC  with  the  following  specihcations: 

•  Length  £  =  8  feet 

•  Tilted  from  vertical  in  the  x  x  2;-plane  by  9t  =  45° 

•  Lifted  over  the  ground  by  he  =  6  ft 

•  Shifted  along  the  x-axis  by  10  ft 

Table  24  lists  the  MiniNEC  description  for  completeness,  where  the  dipole  is  driven 
at  29  MHz.  Figure  40  reports  the  resulting  electric  helds  measured  on  a  hemisphere 
with  a  radius  of  80  ft  or  24  m.  At  29  MHz,  the  wavelength  is  approximately  10  m  so 
that  electric  helds  are  measured  on  a  hemisphere  roughly  2.5  wavelengths  from  the 
antenna. 
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Table  24:  MiniNEC  description  of  the  shifted,  tilted  dipole. 


D :  \antenna\NearToFarField\TiltedDipoleShifted\ShftedTiltedDipole_80  01-28-2013  14:50: 23 

Tilted  Dipole 

45  degrees  elevation  in  the  x-z  plane 
8  long 

6  ft  above  the  ground 
Shifted  10  ft  along  the  x-axis 


GEOMETRY 

Dimensions  in  feet 
Environment:  perfect  ground 


wire  caps  X  Y 

1  none  12.8284  0 

7.1716  0 


Z 

8.8284 

3.1716 


radius 

.001 


segs 


12 


Number  of  wires  =  1 

current  nodes  =  11 


minimum 

Individual  wires 

wire 

value 

segment  length 

1 

.66666 

segment/radius  ratio 

1 

666.66 

radius 

1 

l.E-03 

maximum 
wire  value 

1  .66666 

1  666.66 

1  l.E-03 


ELECTRICAL  DESCRIPTION 
Frequencies  (MHz) 
frequency 

no.  lowest  step 

1  29.  0 


no .  of 
steps 
1 


segment  length  (wavelengths) 
minimum  maximum 

.0196556  .0196556 


Sources 
source  node 
1  6 


sector  magnitude 

1  1. 


phase 

0 


type 

voltage 
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Figure  40:  MiniNEC  simulation  of  the  shifted,  tilted  dipole  at  29  MHz;  electric  helds 
are  measured  at  80  ft. 
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10.1  Synthesis 


Synthesis  requires  extending  the  electric  helds  from  the  upper  hemisphere  to  the 
entire  sphere.  If  77^ (6^,0)  and  77^ (6*,  0)  denote  the  simulated  electric  helds  on  the 
upper  hemisphere,  make  the  even  extension 

Eg{e,(t))  = 

0  <9  <71/2 

7r/2  <9  <77 

and  the  odd  extension 

E^{9,(p)  =  1 

EHOA) 

0  <  0  <  7r/2 

7r/2  <9  <71 

Tables  25,  26,  27,  and  28  report  the  coefficients  estimated  from  Eq  and  sampled 
on  the  upper  hemisphere.  The  hemisphere  is  sampled  using  5°  increments.  There  are 
=  73  samples  in  0  <  0  <  360°  azimuth  and  Ng  =  19  samples  in  0  <  0  <  90° 
zenith.  The  10-ft  offset  of  the  antenna  required  that  the  number  of  terms  increase  to 
N  =  5  over  N  =  3  when  there  is  no  offset. 


Table  25:  \amn\  estimated  from  Eg  and  77^  sampled  on  the  sphere  at  80  ft  for  the 
shifted,  tilted  dipole. 
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Table  26:  |/9mn|  estimated  from  Eg  and  sampled  on  the  sphere  at  80  ft  for  the 
shifted,  tilted  dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

0 

0 

0 

0 

0 

m  = 

1 

2.79  ■  10“® 

2.48  ■  10“® 

1.80  ■  10“® 

6.84-  10“® 

1.20 

■  10-® 

m  = 

2 

5.52  ■  10“® 

2.36  ■  10“® 

2.29-10-^ 

9.26 

■  10“^ 

m  = 

3 

3.75  ■  10-^ 

8.52-10-^ 

5.98 

■  10“® 

m  = 

4 

1.26-  10“® 

4.54 

■  10“® 

m  = 

5 

2.31 

■  10-9 

Table  27:  |7mn|  estimated  from  Eq  and  E^  sampled  on  the  sphere  at  80  ft  for  the 
shifted,  tilted  dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

9.80  ■  10-9 

1.60  ■  10-^ 

5.56  ■  10-® 

1.99-10-^ 

2.51  ■  10-^ 

m  = 

1 

2.34  ■  10-^ 

9.91  ■  10-9 

2.24  ■  10-® 

7.46  ■  10-9 

2.49  ■  10-® 

m  = 

2 

2.79  ■  10-® 

7.16  ■  10-9 

4.55  ■  10-9 

4.60  ■  10-19 

m  = 

3 

5.78  ■  10-9 

2.53  ■  10-1^ 

3.43  ■  10-19 

m  = 

4 

4.58  ■  10-19 

2.42  ■  10-11 

m  = 

5 

9.58  ■  10-11 

Table  28:  \5mn\  estimated  from  Eg  and  E^  sampled  on  the  sphere  at  80  ft  for  the 
shifted,  tilted  dipole. 


n  =  1 

n  =  2 

n  =  3 

n  =  4 

n  =  5 

m  = 

0 

0 

0 

0 

0 

0 

m  = 

1 

9.05  ■  10-® 

3.27-10-^ 

7.71  ■  10-® 

2.27-  10-® 

2.94 

■  10-® 

m  = 

2 

1.07-  10-® 

2.89  ■  10-® 

1.08-  10-® 

1.37 

■  10-® 

m  = 

3 

5.47  ■  10-® 

1.39-  10-® 

5.88 

■10-7 

m  = 

4 

1.57-  10-® 

3.89 

■  10-® 

m  = 

5 

2.93 

■  10-® 
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10.2  Far  Field 


Figure  41  compares  the  MiniNEC  fields  to  Sarkar’s  synthesized  fields  at  a  500-ft  ra¬ 
dius.  The  upper  panels  are  the  sampled  fields  computed  from  MiniNEC.  The  lower 
panels  are  the  synthesized  fields  computed  from  Sarkar’s  expansion  using  the  coeffi¬ 
cients  amm  Pmn,  1mm  ^mn,  estimated  from  the  sampled  fields  on  the  upper  hemisphere 
at  the  80-ft  radius.  These  images  use  the  same  linear  color  scale  with  the  magnitude  of 
the  fields  plotted  in  V/m.  Comparison  between  the  sampled  fields  from  MiniNEC  and 
the  synthesized  field  of  Sarkar  shows  that  differences  in  the  magnitudes  are  relatively 
small. 
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Figure  41:  Comparing  far-field  patterns  of  the  shifted,  tilted  dipole  at  29  MHz.  (upper 
panels)  sampled  fields  computed  from  MiniNEC  at  the  500-ft  radius,  (lower  panels) 
synthesized  fields  computed  at  the  80-ft  radius  and  extended  to  the  500-ft  radius. 
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11 


“V”  Dipole 

Figure  42  is  a  schematic  of  a  half-wavelength,  “V”  dipole  lifted  over  the  ground  plane. 


Figure  42:  “V”  dipole  in  the  y  x  2;-plane  and  lifted  over  a  perfect  ground. 

This  dipole  is  simulated  in  MiniNEC  at  300  MHz: 

•  Each  leg  of  the  “V”  has  length  of  1/4  m 

•  The  angle  between  the  legs  is  Oy  =  45° 

•  The  vertex  is  lifted  over  the  ground  by  he  =  2  m 

Table  29  lists  the  MiniNEC  description  for  completeness,  where  the  dipole  is  driven 
at  300  MHz.  Figure  43  reports  the  resulting  electric  helds  measured  on  a  hemisphere 
with  a  radius  of  80  ft  or  24  m.  At  300  MHz,  the  wavelength  is  approximately  1  m 
so  that  electric  helds  are  measured  on  a  hemisphere  roughly  24  wavelengths  from  the 
antenna. 
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Table  29:  MiniNEC  description  of  the“V”  dipole. 


D :  \antenna\NearToFarField\vDipole\vDipole01_300MHz_80f eet  02-06-2013  07:46: 29 

Vee  dipole:  45  degrees 
2  m  above  perfect  ground 

GEOMETRY 

Dimensions  in  meters 
Environment:  perfect  ground 


wire  caps  X 

Y 

Z 

radius 

segs 

1  none  0 

-.0956709 

2.23097 

.0025 

5 

0 

0 

2. 

2  none  0 

0 

2. 

.0025 

5 

0 

.0956709 

2.23097 

Number  of  wires 

=  2 

current 

nodes  =  9 

minimum 

maximum 

Individual  wires 

wire  value 

wire  value 

segment  length 

1  .05 

1  .05 

segment/radius  rat 

io  1  20. 

1  20. 

radius 

1  2 . 5E- 

03 

1  2.5E-03 

ELECTRICAL  DESCRIPTION 
Frequencies  (MHz) 
frequency 

no.  lowest  step 

1  299.8  0 


no .  of 
steps 
1 


segment  length  (wavelengths) 
minimum  maximum 

.05  .05 


Sources 

source  node  sector  magnitude  phase 

15  11.  0 


type 

voltage 
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Figure  43:  MiniNEC  simulation  of  the  “V”  dipole  at  300  MHz;  electric  fields  are 
measured  at  80  ft. 
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11.1  Synthesis 


Synthesis  requires  extending  the  electric  helds  from  the  upper  hemisphere  to  the 
entire  sphere.  If  Eg{9,(j))  and  E^{9,(1))  denote  the  simulated  electric  helds  on  the 
upper  hemisphere,  make  the  even  extension 


Ee{9,(^) 

and  the  odd  extension 

E^i9,(p)  = 


0  <  0  <  7r/2 

7r/2  <  9  <  TT 

El(e,4>) 

0  <  0  <  7r/2 

1 

1 

7r/2  <  9  <  71 

Figure  44  reports  Sarkar’s  coefficients  estimated  from  Eq  and  E^  sampled  on  upper 
hemisphere  of  an  80-ft  sphere.  This  hemisphere  is  sampled  using  5°  increments.  There 
are  =  73  samples  in  azimuth  (0  <  0  <  360°)  and  Nq  =  19  samples  in  zenith 
(0  <  0  <  90°). 


2  4  6  8  10  12  14  16  18 


2  4  6  8  10  12  14  16  18 


2  4  6  8  10  12  14  16  18 


Figure  44:  Sarkar’s  coefficients  estimated  from  Eq  and  E^  sampled  on  the  sphere  at 
80  ft  for  the  “V”  dipole. 

The  coefficients  are  plotted  in  dB  (i.e.,  10  x  log;^Q(|a(m,  n)  |))  so  that  each  decade 
corresponds  to  a  power  of  10.  The  plot  shows  contributions  from  (3  and  7.  The  slow 
“roll-off”  of  (3  implies  that  convergence  will  be  slow. 
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For  example,  Figure  45  reports  on  using  iV  =  11  coefficients  to  synthesize  the 
electric  helds  at  80  ft.  The  upper  panels  are  the  electric  helds  computed  from  Mini- 
NEC  at  a  radius  of  80  ft  and  are  the  “ground  truth.”  The  lower  panels  are  the 
synthesized  electric  helds  at  80  ft.  The  synthesis  starts  by  sampling  both  helds  on 
the  80-ft  upper  hemisphere  at  5°  increments,  extending  both  helds  to  the  whole  80- 
ft  sphere,  estimating  Sarkar’s  coefficients  from  the  samples,  and  hnally,  using  the 
estimated  coefficients  to  synthesize  the  held  at  r  =  80  ft.  The  lower  panels  show  that 
the  gross  structure  of  the  helds  are  approximated  while  the  magnitudes  needs  more 
terms  to  “hll-in.” 
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Figure  45:  Near-held  synthesis  using  N  =  11,  M  = 


0.6 

0.4  .E 

0.2 

0 


3. 


Figures  46  and  47  show  that  increasing  the  number  of  terms  produces  synthetic 
helds  with  increasing  accuracy.  The  hnal  plot  is  considered  “good  enough”  to  extend 
to  the  far  held. 
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Figure  46:  Near-field  synthesis  using  TV  =  15,  M  = 


0.6 
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Figure  47:  Near-field  synthesis  using  =  19,  M  =  3. 
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11.2  Far  Field 


Figure  48  reports  on  the  far-held  synthesis  at  500  ft.  The  upper  panels  are  the  electric 
helds  computed  from  Mini-NEC  at  500  ft.  The  lower  panels  are  the  corresponding 
synthetic  held  at  500  ft.  These  synthetic  helds  were  produced  by  sampling  the  Mini- 
NEC  helds  at  r  =  80  ft  on  upper  hemisphere  at  5°  increments,  extending  both  helds 
to  the  whole  80-ft  sphere,  estimating  Sarkar’s  coefficients  from  these  samples,  and 
hnally,  using  these  estimated  coefficients  to  synthesize  the  held  at  r  =  500  ft.  The 
lower  panels  show  that  the  =  19  terms  produces  an  approximation  that  can  support 
a  gain  computations. 
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Figure  48:  “V”  dipole  synthesis  at  r  =  500  ft  from  Eq  and  sampled  on  the  upper 
hemisphere  at  80  ft. 
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11.3  Far-Field  Gains 


The  radiation  intensity  is  defined  for  the  far  field  as  [4,  Equation  2-12a] 

2 

U(e,<i>)  :=  ^  {|£»(«,^)|" +!£♦(», ^)iq  [Watts], 

Figure  49  compares  the  radiation  intensity  for  the  “V”  dipole  between  the  Mini- 
NEC  “ground  truth”  (upper  panel)  and  the  a  synthetic  approximation  constructed 
from  =  19  terms  (lower  panel).  The  range  for  both  fields  is  500  ft.  At  300 
MHz,  this  corresponds  to  152  wavelengths.  The  synthetic  radiation  intensity  exhibits 
discrepancies  for  increasing  zenith  or  decreasing  elevation  angle.  For  zenith  of  0°,  the 
synthesis  algorithms  used  in  this  report  set  the  field  to  zero 


Figure  49:  “V”  dipole  radiation  intensity  at  r  =  500  ft  from  Eg  and  sampled  on 
the  upper  hemisphere  at  80  ft. 
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12  Phase  2  Tasking 


Phase  1  assessed  an  implementation  of  Sarkar’s  near-to-far  transformation  [12]  on 
analytic  antenna  models  and  Mini-NEC  antenna  simulations.  For  both  classes  of  an¬ 
tenna  models,  Phase  1  verihed  the  implementation  recovers  the  far  held  from  the  near 
held — provided  the  spherical  expansion  includes  a  sufficient  number  of  terms.  That  is, 
the  truncation  error  governs  the  near-to-far  transform  rather  than  the  discretization 
error — assuming  noise-free  antenna  models  (see  the  discussion  for  Equation  2). 

Phase  2  applies  Sarkar’s  near-to-far  transformation  to  measurements  taken  at  SSC 
Pacihc’s  Antenna  Pattern  Range.  These  electric  held  measurements  are  degraded  by 
additive  noise,  in-band  RE  interference,  and  the  mechanical  and  electrical  variations 
in  the  Antenna  Pattern  Range.  Therefore,  an  initial  task  for  Phase  2  is  to  assess  the 
robustness  of  the  near-to-far  transformation  with  respect  to  the  noise-corrupted  mea¬ 
surements.  Figure  50  shows  the  candidate  reference  antenna  placed  on  the  Antenna 
Pattern  Range.  This  antenna  consists  of  a  monopole  and  a  loop  above  the  ground 
plane. 


Figure  50:  Reference  antenna  in  the  SSC  Pacihc  Antenna  Pattern  Range. 

Because  this  antenna  is  amenable  to  modeling,  comparison  of  the  ideal  model 
against  the  measured  antenna  is  possible  and  will  produce  a  noise  model  for  the 
spherical  coefficients.  Although  the  near-to-far  transformation  is  not  expected  to 
amplify  the  noise  in  the  far  held,  the  accuracy  of  the  spherical  coefficients  may  degrade 
the  asymptotic  estimates  of  the  truncation  error.  The  noise  model  can  set  a  truncation 
limit — increasing  the  number  terms  in  the  expansion  may  not  improve  the  accuracy 
because  the  higher-order  coefficients  are  lost  in  the  noise. 
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Classical  statistical  smoothing  is  a  standard  approach  to  reduce  additive  noise  but 
will  not  mitigate  system  measurement  errors.  Assessing  the  payoff  for  this  smoothing 
is  a  necessary  Phase  2  task.  A  more  promising  approach  uses  model  identihcation  to 
mitigate  noise.  In  this  approach,  an  antenna  model  is  htted  to  the  noisy  measured  near 
held.  For  example,  measurements  of  the  candidate  monopole  may  be  htted  by  varying 
antenna  length,  radius,  orientation,  and  material.  Once  a  model  is  determined  that 
best  “explains”  the  measurements,  this  model  is  then  used  to  compute  a  “noise-free” 
far  held  and  extrapolate  to  other  frequencies. 
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